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ABSTRACT 

We present the results of a high-resolution spectroscopic monitoring program of the eclipsing 
pre-main-sequence star KH 15D that reveal it to be a single- line spectroscopic binary. We find that 
the best-fit Keplerian model has a period P — 48.38 days, which is nearly identical to the photometric 
period. Thus, we find the best explanation for the periodic dimming of KH 15D is that the binary 
motion carries the currently visible star alternately above and below the edge of an obscuring 
cloud. The data are consistent with the models involving an inclined circumstellar disk, as recently 
proposed by Winn et al. (2004) and Chiang & Murray-Clay (2004). We show that the mass ratio 
expected from models of PMS evolution, together with the mass constraints for the visible star, re- 
strict the orbital eccentricity to 0.68 < e < 0.80 and the mass function to 0.125 < Fm / sin 3 i < 0.5 M . 



Subject Headings: stars: pre-main-sequence — stars: individual (KH 15D) 
troscopic — techniques: radial velocities — circumstellar matter 
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binaries: spec- 



1. INTRODUCTION 

KH 15D is a K6-K7 pre-main-sequence star 
that exhibits drama tic photometric variability 
l)Kearns fc Herbstl 119981) . Every 48 days, the star's 
brightness dims by 3.5 magnitudes and remains in this 
faint state for nearly half of the photometric period. 
These deep brightness minima are acco mpanied by a 
slight blueing of the star's color indices l)Herbst et al.l 
120021 hereafter He02), littl e or no change in spectral 
type IHamilt on et al.l 120011) a nd an increase in linear 
polarization ijAgol et al.l 12001 . This implies that the 
star is completely eclipsed by an optically thick, ex- 
tended collection of dust grains, possibly in the form of 
a circumstellar disk. If this is the case, the serendipitous 
alignment of the KH 15D star/disk system may provide 
insights into the evolution of young stars and their 
interactions with their circumstellar environments. 

While KH 15D has periodic eclipses, it cannot be an 
ordinary eclipsing binary because of the long duration of 
minimum light. Recent theories postulate that there is, 
nonetheless, a currently unseen binary companion to the 
visible star. In a study of archival photographic plates, 
Uohnson fc Wind 1)20041 hereafter JW04) discovered that 
the historical light curve of KH 15D is similar to the mod- 
ern light curve but appears to be diluted b y light from 
a seco nd star. Motivated by these findings. iWinn et al.l 
l)2004l hereafter W04) constructed a model composed of 

1 Based on observations obtained at the W. M. Keck Observa- 
tory, which is operated jointly by the University of California and 
the California Institute of Technology; Las Campanas Observatory 
of the Carnegie Institution with the Magellan II Clay telescope ; 
and McDonald Observatory of the University of Texas at Austin. 

2 Department of Astronomy, University of California, Mail Code 
3411, Berkeley, CA 94720 

3 Visiting astronomer, McDonald Observatory which is operated 
by The University of Texas at Austin. 

4 Department of Physics and Astronomy, Rice University, Hous- 
ton, TX 77005. 

5 Van Vleck Observatory, Wesleyan University, Middletown, CT 
06459 



a binary system with the orbital plane inclined with re- 
spect to the edge of an optically thick screen. As the 
two stars orbit one another, the reflex motion carries one 
star alternately above and below the edge of an opaque 
screen, causing the eclipses. The long-term evolution of 
the light curve is reproduced by allowing the screen to 
move slowly across the binary orbit — a feature of the 
model that led the authors to envision the screen as a 
processi ng, circumbinary disk. A sim ilar model is pro- 
posed bv lChiang fc Murrav-Clavl 1)2004 hereafter CM04) 
who also envision the opaque screen as an inclined, pro- 
cessing circumbinary disk, or "ring," with an inner edge 
truncated by tidal interactions with the binary and an 
outer edge possibly shepherded by a planetary compan- 
ion. Both models provide explanations of the unique 
features of KH 15D's light curve. The W04 model makes 
quantitative predictions about the orbital parameters of 
the binary system, while the CM04 model provides a 
physical description of the circumbinary ring. 

A fundamental question that has not yet been an- 
swered is whether KH 1 5D is a single or multi ple stellar 
system. He02 (see also IHamilton et alJ 120031 hereafter 
Ha03) first searched for evidence of orbital companions 
using high-resolution VLT spectra and reported a radial 
velocity change of +3.3 ± 0.6 km s _1 over two widely 
spaced epochs. However, one of the measurements was 
made during egress when there was a strong possibility of 
contamination from scattered light or line profile distor- 
tions from a sharp-edged cloud only partially occulting 
the stellar disk. 

Over the past two years we have conducted a high- 
resolution, multi-site spectroscopic monitoring campaign 
to determine whether or not KH 15D exhibits orbital mo- 
tion indicative of a multiple system. Here we present the 
results of our study, which show that KH 15D undergoes 
significant radial velocity variations. The variations are 
consistent with a binary companion with an orbital pe- 
riod equal to the 48-day photometric period. In we 
summarize our observations and reduction procedures. 
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TABLE 1 

Spectroscopic Observations of KH 15D 



UT Date 


Telescope/ 
Instrument 


AA [A] 


R 


J.D.-2.4 x 10 6 


v r [km/s] 


9^ Drt 9009 




^qoo fi9nn 


70 000 


^9^79 ^74 


1.7 rt 0.2 a 


28 Oct 2002 


Keck/HIRES 


3900-6200 


70,000 


52575.511 


3.0 ± 0.3 a 


16 Dec 2002 


Magellan/MIKE 


4500-6800 


45,000 


52624.738 


3.1 ± 0.2 b 


14 Jan 2003 


Keck/HIRES 


3900-6200 


70,000 


52653.468 


3.3 ± 0.2 a 


R Feh 900^ 






70 000 


^9fi78 400 


9.0 i 0.5 b 


9 Feb 2003 


Keck/HIRES 


4200-6600 


70 000 


52679 41 


11.5 ± 0.5 b 


9 Mar 2003 


Magellan/MIKE 


4900-8700 


22,000 


52707.511 


0.8±0.2 b 


3 Nov 2003 


Keck/HIRES 


4200-6600 


30,000 


52946.584 


1.2±0.4 b 


4 Nov 2003 


Keck/HIRES 


4200-6600 


30,000 


52947.594 


1.8±0.3 b 


4 Jan 2004 


McDonald/CE 


5600-6900 


30,000 


53008.762 


1.4±0.6 b 


5 Jan 2004 


McDonald/CE 


5600-6900 


30,000 


53009.852 


1.7±0.4 b 


10 Jan 2004 


Keck/HIRES 


3900-6200 


70,000 


53014.320 


5.4±0.3 a 


10 Jan 2004 


McDonald/CE 


5600-6900 


30,000 


53014.760 


5.8±0.4 b 


11 Jan 2004 


Keck/HIRES 


3900-6200 


70,000 


53015.506 


7.0±0.2 b 


9 Feb 2004 


Keck/HIRES 


3900-6200 


70,000 


53044.834 


1.8±0.2 a 


10 Feb 2004 


Keck/HIRES 


3900-6200 


70,000 


53045.828 


1.3±0.2 a 



a Multiple reference stars were used to calculate the radial velocity. The reported velocity is the mean velocity measured from the various 
target-reference star pairs. The estimated uncertainty is the standard deviation of the mean velocity. See Tablc|5]for a listing of individual 
reference stars and radial velocity measurements for each target-reference star pair. 

b Only one reference star was used in the radial velocity measurement. The reported velocity is the mean velocity from all echelle orders. 
The estimated uncertainty is the standard deviation of the mean velocity. 



The radial velocity measurements and best-fit Keplerian 
orbital parameters are presented in fJ3J We conclude in 
2] with a discussion of our findings and implications for 
existing models of the KH 15D eclipse mechanism. 

2. DATA 

During the 2002-2003 and 2003-2004 observing sea- 
sons we made 16 observations of KH 15D at maximum 
light including 11 spectra with the Keck 10- meter tele- 
scope and HIRES echelle spectrometer; 2 spectra with 
the 6.5-meter Magellan II (Clay) telescope and the MIKE 
echelle spectrometer; and 3 spectra with the 2.1-meter 
Otto Struve Telescope at McDonald Observatory with 
the Sandiford Cassegrain Echelle Spectrometer (CE). 
The observations are summarized in Table ^ Addi- 
tional Keck/HIRES spectra were obtained during min- 
imum light. However, since it is unlikely that the star's 
photosphere is visible through the obscuring material, we 
decided to exclude from our analysis spectra that were 
obtained within 10 days of mid-eclipse. In addition to 
our own measurements, we also include the out-of-eclipse 
radial velocity measurement, v r = 9.0 ± 0.2 km s , re- 
ported by (Ha03) based on their VLT/UVES spectra. 

Most of the Keck/HIRES observations were made as 
part of the California & Carnegie Planet Search. 6 For 
these observations, the relatively faint apparent magni- 
tude of KH 15D (V — 16) at maximum light precluded 
the the use of the iodine cell to establish a wavelength 
scale. The cell was therefore removed from the light path 
during observations of KH 15D to increase the through- 
put of the spectrometer. The raw CCD frames from all 
telescopes were reduced using reduction packages writ- 
ten in IDL. The details of the reduction procedures are 
fundamentally identical to the algorithm presented by 
iValentil l)1994|) . After bias-subtraction, each echelle frame 
is divided by a normalized median flat-field image. Or- 
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der definition is performed using a bright star or flat-field 
exposure, and scattered light is removed by fitting a two- 
dimensional B-spline to the inter-order regions and inter- 
polating across each spectral order. After the scattered- 
light is subtracted, each order is rectified, sky-subtracted 
and summed in the cross-dispersion direction to form the 
final one-dimensional spectrum. 

Instead of a summation in the cross-dispersion direc- 
tion, the rectified orders of the McDonald spectra are 
reduced to one-dimensional spec tra using the optima l ex- 
traction algorithm described bv lHinkle et al.l (|2000^l . In 
the case of the MIKE reductions, the standard code is 
modified to correct for the tilt of the spectrometer en- 
trance slit with respect to the CCD columns. The cor- 
rection of the slit tilt is necessary because the sky sub- 
traction algorithm we employ requires that the projected 
slit image lie parallel to the detector columns. 

The radial velocity of KH 15D relative to the Solar 
System barycenter is measured from each spectral ob- 
servation by means of a cross-correlation analysis. For 
the spectra of KH 15D obtained as part of the Planet 
Search observing program, the program stars observed 
on each night provided an extensive selection of refer- 
en ce stars with known b arycentric radial velocities listed 
in iNidever et al.l l)2002|) . We selected reference stars 
that were observed within 30 minutes of KH 15D and 
with spectral types ranging from M0 to G5. Since the 
Planet Search target stars are observed through an io- 
dine cell, orders containing iodine absorption lines are 
avoided in the analysis. For KH 15D observations ob- 
tained as part of programs other than the Planet Search, 
a single K-type reference star observed either before or 
after KH 15D is used as the reference star. 

The cross-correlation of each KH 15D spectrum with 
respect to the reference spectra is performed using cus- 
tom routines written in IDL. The procedure involves first 
rcbinning each one-dimensional spectral order onto a new 
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TABLE 2 

Absolute Radial Velocity Measurements of KH 15D From 



Each Reference Star 


J.D.-2.4 x 10 6 


Reference 


v r [km/s] a 




Star HD# 


rar TO If f? A 

52572. 5 74 


A Afl7fl 

4U979 


1.6 ±0.1 




56274 


2.1 ±0.1 




52265 


n 1 I n no 

2.1 ± 0.08 




49674 


"i i I n nn 

1.3 ± 0.09 




66171 


1.2 i 0.2 


52575.511 


31966 


2.9 1 0.2 




36395 


a r\ _l_ n no 
4.0 ± 0.08 




47157 


3.2 i 0.2 




37124 


2.2 ± 0.1 




a nn7n 

4uy (» 


2.7 ± 0.3 


rr O/" 1 O A TO o 


10700 


3.1 ± 0.2 


020DO.40O 


31966 


3.6 ±0.1 




42618 


2.7 ±0.1 




47157 


Q 1 _L n (17 

3.1 ± 0.0/ 




99109 


q a _l_ n nn 
3.0 ± 0.U9 




nn /i no 

99492 


3.4 ±0.1 


roc7o /inn 


obUUb 


9.0 ± 0.5 




OOUUo 


11.5 ± 0.5 


52 (U f .511 


ooUUo 


0.8 ± 0.2 


5zy4o.oo4 


36006 


1.2 ± 0.4 


c on d7 r* n -i 

52947.594 


36006 


1.8 ± 0.3 


tronno 7ro 

53008.762 


36006 


1.4 ± 0.6 


53009.852 


17660 


1.7 ± 0.4 


53014.320 


45350 


5.1 ±0.2 




18830 


5.6 ±0.3 


53014.760 


17660 


5.8 ±0.4 


53015.506 


45350 


7.0 ±0.2 


53044.834 


50499 


2.0 ±0.2 




46375 


1.7 ±0.1 




50281 


1.8 ±0.2 




46375 


1.9 ±0.2 


53045.828 


50499 


1.6 ±0.1 




76909 


1.3 ±0.1 




73667 


0.9 ±0.1 




72673 


1.2 ±0.1 



a Velocities are the mean value from all echelle orders used in the 
cross-correlation analysis. Estimated uncertainties are the stan- 
dard deviation of the mean velocity from all echelle orders. 



wavelength scale that is linear in log A. This ensures that 
each pixel in the rebinned spectrum represents a veloc- 
ity interval that is un iform over the entire spectral order 
l)Tonrv fc Davisiri979fl . Regions containing telluric lines, 
strong emission features, and CCD defects are masked 
out and each spectral order is cross-correlated with re- 
spect to the corresponding order of the reference spec- 
trum. Each spectral order thus yields an independent 
measurement of the radial velocity of KH 15D relative 
to the reference star. The average of the ensemble set of 
velocities from all orders is then adopted as the relative 
radial velocity of KH 15D for a given epoch. 

The relative radial velocities from each night are con- 
verted into absolute barycentric radial velocities using 
the relation 



Av + (BC kh - BC re f) + V re f. 



(1) 



In Eqn. JQ), Av is the relative velocity from the cross- 
correlation analysis; BCkh and BC re f are the barycentric 
corrections for KH 15D and reference star, respectively; 
and v re f is the absolute ba rycentric radia l veloc ity of 
the reference star as listed in lNidever et alJ l)2002f) . The 
mean measured velocities from each epoch are listed in 
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Fig. 1. — The radial velocity of KH 15D as a function of time. 
The dashed line is the best-fit Keplerian orbit with the eccentricity 
fixed at e = 0.74 — the mean value allowed by our orbit constraints 
(see 



the final column of Table Table lists the velocities 
obtained from our cross-correlation analysis for each in- 
dividual reference star on each night. 

For nights when only one reference spectrum is avail- 
able, the uncertainty is estimated using the standard de- 
viation of the mean radial velocity measured from all 
echelle orders. For nights with multiple reference spectra, 
the standard deviation of the velocities computed from 
each target-reference pair is adopted as the uncertainty. 
This latter method of error estimation typically yields a 
larger value than the order-to-order scatter seen for an 
individual observation (0.2 < a v < 0.6 km s _1 compared 
to 0.1 < a v < 0.3). The larger scatter in velocity mea- 
surements among the target-reference star combinations 
is likely due to external systematic effects such as small 
shifts in the position of the CCD throughout the night 
caused by thermal relaxation of the detector mount. 
Thus, the standard deviation measured from the echelle 
orders in a single observation likely underestimates the 
true uncertainty by approximately a factor of two. 

3. RESULTS 

Figure ^ illustrates how the radial velocity of 
KH 15D varies temporally out of eclipse over a range 
of 10.7 km s _1 . These data indicate that there must 
be an unseen star in the system, as had previously been 
inferred from the analysis of the historical light curve 
(JW04). Hereafter we will adopt the naming convention 
of W04 and refer to the currently visible star as A and 
the hidden companion as B. 

3.1. Orbit Solution 

Using a nonlinear least-squares algorithm, we found 
a best-fit model orbit with a period P — 48.38 days, 
eccentricity e > 0.27 and velocity semiamplitude K > 
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Fig. 2. — The radial velocity of KH 15D as a function of orbital phase. The solid line is the best-fit Keplerian orbit with the eccentricity 
fixed at e = 0.74 — the mean value allowed by our orbit constraints (see ilij.21 . The fit has reduced chi-squared \/x£ = 1.5. The vertical 
dotted lines represent the approximate phases of ingress (left) and egress (right) based on the He02 ephemeris. The bottom panel shows 
the 2001-2002 /-band light curve. The phased photometric measurements have been binned for clarity using phase intervals of Acj> = 0.015. 



8.15 km s _1 . The other orbit parameters are listed in 
Table 03 Figure shows a plot of radial velocity versus 
orbital phase for P = 48.38 days. The rms scatter of the 
fit residuals is 0.38 km s _1 and the reduced \/xJ = 1.5. 
The vertical lines at <f> = 0.33 and (f> = 0.73 denote the 
approximate phases of ingress and egress, respectively, 
based on the He02 ephemeris. We have no radial veloc- 
ity measurements between these phases because of the 
eclipse of star A. 

We find that the orbit solution is not well constrained 
due to the lack of data near periapse, which allows the 
velocity semi-amplitude of the orbit solution to compen- 
sate for changes in the eccentricity. As such, we were 
able to obtain reasonable fits by fixing the eccentricity 
at values e > 0.27, with each solution yielding different 
values of the orbit parameters and values of \fxl that 
are equivalent at the 97.5% confidence level (based on 
the 9 degrees of freedom in the fit). 

Figure|21shows how and the orbit parameters vary 
as a function of eccentricity. While the period varies little 
over a wide range of eccentricities, it is apparent that the 
radial velocity data alone provide poor constraints for 
the other orbit parameters. However, the implied mass 
ratio from other KH 15D observations, together with the 
measured mass limits of star A, can be used to place 



limits on the orbital eccentricity, as we now show. 

3.2. Orbit Parameter Constraints 

JW04 determined from photometric measurements 
of archival plates that the out-of-eclipse magnitude of 
KH 15D was 0.9 mag brighter (at /-band) 40 years ago 
compared to the modern bright state. The two most 
probable explanations for the brighter apparent magni- 
tude in the past are that either both stars were visible, 
or that star B alone was visible. If both stars were vis- 
ible, then Lb/ La = 1.3. If only B were visible, then 
Lb I La — 2.3. In either case, the condition L b /La > 1 
must hold. 

Multi-color photometric measurements obtained dur- 
ing minimum light show a slight blueing of the color in- 
dices (Ha 03) compared to t he colors at maximum light. 
Similarly, Ago l et al.1 l)2004|) measure slightly bluer col- 
ors during eclipse from their low-resolution spectropo- 
larimetric observations compared to their out-of-eclipse 
observations. One possibility for the bluer color indices 
during eclipse is that the scattered light is dominated by 
Rayleigh scattering. However, He02 show that there is 
no reddening of the light from star A during ingress and 
egress. This suggests that the opacity of the occulting 
material is wavelength-independent and that the bluer 
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0.10 0.20 0.30 



0.40 0.50 0.60 
Orbital Eccentricity 



0.70 0.80 



0.90 



Fig. 3. — The dependence of the best-fit orbit parameters (top 
5 panels) and reduced chi-squared (lower panel) on the the orbital 
eccentricity. The ordinate values are generated by fixing the orbital 
eccentricity in the least-squares fit to the radial velocity data. The 
vertical dashed line indicates the lower limits placed on the fit 
parameters based on the 97.5% confidence level for y/ xt with 9 
degrees of freedom. The solid vertical lines denote t he up per and 
lower limits placed on the eccentricity as discussed in it3.2l (see also 
Figure|lJ, and the fits generated for e = 0.74 are denoted by solid 
dots with error bars. 



M A (e,M A /M B ) = ¥± 
Mb 



1 



Ma 
Mb 



F M (e,K,P) 



sin 3 i 



(3) 

The mass function Fm/ sin 3 i in Eqn. (j3Jl is calculated 
from the parameters of the best-fit orbit to the radial 
velocity data with the value of e fixed in the fitting 
procedure. Figure 0] shows plots of Ma versus e for 
Ma/Mb = 0.7, 0.8, 0.9 and 1.0. Also shown is the 
value of Ma/Mb — 1-6 obtained by W04 based on their 
model fit to the KH 15D light curve (see ©. Since the 
mass of the visible sta r is known to fall within the limits 
0.5 < M A < 1.0 M Q l|Flaccomio et al.lll999t iPark et alJ 
2000) , it can be seen by inspection of Figure 0] that ac- 
ceptable choices of the eccentricity fall within the range 
0.68 < e < 0.80. 

The resulting range of acceptable orbit parameter val- 
ues is shown in Figure |3] between the solid, vertical lines 
in each panel. Column two of TableOllists the best-fit pa- 
rameter values obtained using the mean acceptable value 
of the eccentricity, e = 0.74 (see also Figure|5J). The un- 
certainties in the fit parameters were estimated using a 
Monte Carlo simulation. We generated 10 3 statistical re- 
alizations of the velocity time series assuming the errors 
are normally distributed with standard deviations equal 
to the measurement uncertainties. The estimated uncer- 
tainties are also displayed in Figure as error bars on 
the points located at e = 0.74. 

3.3. Periodicity 

As discussed in ti3.ll we find the best-fit Keplerian has 
an orbital period of P = 48.38 ± 0.01 days. He02 report 
a period of P — 48.35 ±0.02 days from their photometric 
monitoring and JW04 report P = 48.42 ±0.02 days from 
a periodogram analysis of archival photographic plates 
photometry. Thus, our measured orbital period differs 



colors are due to a bluer object. Therefore the temper- 
ature of B must be hotter than the temperature of A, 
assuming both stars contribute nearly equally to the scat- 
tered component of the light observed during the eclipse 
of A. However, since the colors are only ^0.1 mag bluer 
at minimum light (He02), Tb cannot be much larger than 
T A . 

For most low- mass (Af» < 1 .0 M^) pre-main- 
sequence evolutionary models (e.g iChabrier &: Baraffd 
119971 iD'Antona fc MazziteTIH Il997|i the stellar mass is 
monotonic with both luminosity and temperature. For 
coeval stars on their Hayashi tracks, Tb > T A and 
L b /La > 1 imply that Mb > Ma or, in terms of the 
mass ratio, Ma/Mb < 1. 

The mass function of a Keplerian orbit can be ex- 
pressed as 




0.4 0.5 0.6 0.7 

eccentricity 



0.8 



F M (e,K,P) = 



K 3 P(l-e 2 )ism 3 i 
2ttG 



M§ sin 3 i 
(Ma + Mb) 



where P is the orbital period and K is the velocity semi- 
amplitude. Solving Eqn. J5J for Ma yields 



Fig. 4. — The mass of the visible component of KH 15D, 
Mj\, versus orbital eccentricity for five choices of the mass ratio 
M A /M B . The dotted lines show for 0.5 < M A < 1.0 M Q and 
Ma/Mb ^ 1-0, acceptable values of the eccentricity fall within the 
range 0.68 < e < 0.80. The inclination is assumed to be i = 90° 
(sin^i = 0). 
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TABLE 3 
Optimized Model Parameters 



Parameter 


Acceptable Fit to 


Acceptable Fit With 


W04 




Velocities' 1 


e = 0.74 


Model 1 



P [day] 48.38 < P < 48.40 48.38(0.01) 48.35 (fixed) 

e > 0.27 0.74 (fixed) 0.7 

w [deg] -20 < ui < 1 -2 (1) -7.2 

7 b [km s" 1 ] +6.9 < 7 < +22.5 +14.7 (0.9) +15.5 

asini [AU] < 0.29 0.21 (0.5) 0.18 

K [km s" 1 ] > 8 53 (4) 27.5 

F AI /sin 3 i 2.4 X 10- 3 < F AI /sin 3 i < 2.2 0.24 (0.05) 3.8 X 10~ 2 

T p c [J.D.] 2452251.0 (0.6) 

Fit rms [km s" 1 ] 0.38 0.38 

Reduced \J_ x„ 1.5 1.5 - 

a The limits on the fit parameters correspond to the 97.5% confidence upper limit on v^X? ( sce Figure[3J and the condition e < 1. 
b Radial velocity of the binary C.O.M. with respect to the Solar System barycenter 
c Time of periastron passage 



by 1.5cr and 2a from the modern and historical photo- 
metric periods, respectively. Since the evolving shape 
of the light curve may complicate the accurate determi- 
nation of the photometric period, we conclude that the 
orbital period and photometric period of KH 15D are 
indistinguishable within measurement errors. 

Hc02 noted that the shape of the light curve varies 
from eclipse to eclipse in such a way as to suggest a 96- 
day period for the system, with each cycle containing 
two brightness minima. In order to investigate this al- 
ternative periodicity, we fixed the period in our fit at 
P = 96.76 days and obtained a reduced -\/xS of 47.5. 
Based on this poor quality of fit, we find that a 96-day 
period is incompatible with our radial velocity measure- 
ments. 

4. DISCUSSION 

The results of our spectroscopic monitoring campaign 
show that KH 15D is a spectroscopic binary system. The 
observed radial velocities are consonant with a stellar 
companion with an orbital period equal to the photo- 
metric period. We now discuss the implications of the 
binarity of KH 15D for current models of the photomet- 
ric variability mechanism. 

Existing models of the KH 15D eclipse mechanism fall 
into two classes differentiated by whether it is the orbital 
motion of the star or a feature of a circumstellar disk that 
causes the photometric variability. The first class of mod- 
els posits the existence of a single star surrounded by a 
circumstellar disk containing a nonaxisymmetric density 
enhancement or alternatively, a warp. As the disk fea- 
ture orbits the star with a 48-day period, it periodically 
blocks the line-of-sight to the stellar surface. Based on 
their spectropolarimetric observations that sh owed an in- 
crease in polarization during minimum light. lAgol et all 
l)2004[) developed a model of a warped disk with an ex- 
tended atmosphere and obtained a reasonable fit to the 
20 01-2002 light curve. A similar analysis was performed 
bv lBarge fc Vito3 l|2003D using a model involving a large 
dusty vortex. 

While these single-star models are able to match the 
light curve of KH 15D at one point in time, they do not 
adequately explain the temporal changes in the observed 
light curve. Perhaps more importantly, it is not clear if 



a binary companion to KH 15D is compatible with an 
eclipsing disk feature. If a density enhancement in a cir- 
cumstellar disk is responsible for the 48-day photometric 
period, the feature must orbit the central star at a dis- 
tance of ~ 0.22 AU, assuming a Keplerian orbit. This 
geometry would restrict the stellar companion to orbit 
either at a distance less than 0.22 AU or beyond the 
outer extent of the circumstellar disk. However, since 
our best-fit orbit yields a period of 48.38 ±0.01 days and 
a sin i = 0.21 AU (assuming e = 0.74) for the binary com- 
panion, neither of these scenarios seems plausible. Thus, 
it does not seem possible for both an eclipsing disk fea- 
ture and a stellar-mass companion to coexist with the 
same orbital period. 

In the second class of models, proposed by W04 and 
CM04, there exists an unseen binary companion of com- 
parable mass to the visible K-type star seen today. Sur- 
rounding the two stars is a circumbinary disk viewed 
nearly edge-on. The orbital plane of the two stars is 
tilted at a small angle with respect to the disk plane and 
the present day light curve is produced as the reflex mo- 
tion of Star A carries it above and below the disk plane. 

Independent evidence for the binary nature of 
KH 15D has emerged from studies of archival photo- 
graphic plates. In a study of photographic plates ob- 
tained from Asiago Observatory, JW04 show that the 
apparent magnitude of KH 15D was variable from 1968 
to 1983, but the light curve from this epoch was markedly 
different from the one observed today. The bright state 
was nearly a factor of 2 brighter in the past and the 
eclipse depth was a factor of 5 shallower. These find- 
ings can be explained by invoking the presence of a sec- 
ond star that was visible in the past but is unseen today 
(JW04). 

Motivated by these findings, W04 constructed a model 
of an eccentric binary with a fraction of the orbital plane 
obscured by an opaque screen (or circumbinary disk) and 
found a quantitative orbital solution by fitting simultane- 
ously to the 2001-2002 photometry of He02, the historical 
photometry of JW04 and the radial velocity measure- 
ments of Ha03. Based on a preliminary investigation of 
the radial velocities presented here, CM04 independently 
used physical arguments to propose a similar model of 
an eccentric binary surrounded by a nodally precess- 
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ing circumbinary ring. As the ring precesses, the light 
curve gradually changes from the one recovered from the 
archival plates to the shape seen today. Perhaps the 
greatest advantage of the two-star models is their ability 
to explain not only the present-day light curve, but also 
its evolution over the past half century. 

In addition to explaining the photometric phenomenol- 
ogy of KH 15D, the two-star models make predictions 
about the nature of the binary orbit. Both CM04 and 
W04 predict that (1) periastron passage occurs during 
minimum light, (2) the orbital companion has a mass 
comparable to the currently visible star, and (3) the bi- 
nary orbit is highly eccentric. These predictions are pre- 
cisely what we find from our orbit solution. 

The W04 model makes additional, quantitative predic- 
tions about the orbital parameters of the binary. For an 
assumed fixed period of 48.35 days, the model produces 
a velocity semi-amplitude of 27.5 km s -1 , eccentricity 
e = 0.7, a mass ratio Ma/Mb = 1.6, argument of peri- 
center oj = —7°. 2 and a center of mass radial velocity 
7 = +15.5 km s _1 . As can be seen in TableOl the W04 
Model 1 predictions agree well with our orbital solution 
with fixed eccentricity, e = 0.74 (see jjHUjl - We note here 
that our center of mass velocity for the binary system of 
7 = +14.7 ± 0.9km s" 1 clearly rules out W04 Model 2, 
which predicts 7 = +5.7 kms and negative velocity at 
periapse. 

The W04 model also predicts the inclination of the 
binary orbit, a property of the system that our radial 
velocities do not. W04 Model 1 predicts i = 84°. 6 or 
sini = 0.996. Similarly, using the geometry of the cir- 
cumbinary ring proposed by CM04, an upper limit on 
the inclination of i < 80° can be assigned to the binary 
orbit based on the time lag between periastron passage 
and mid-eclipse (Chiang, private communication). 

We note that the W04 model has the peculiar feature 
that the less massive star is the more luminous star. Us- 
ing Equation^ the values of P, K and e produced by the 
W04 model lead to a mass function Fm — 0.038 M© and 
Ma = 0.41 Mq. This mass is sig nificantly less than the 
lower limit of 0.5 M measured bvlPark et alJ (I2 000fl an d 
the value of 0.6 M© measured by Flac comio et al.l fl999). 
Figure|Ushows a plot of Ma as a function of e for the W04 
mass ratio Mb /Ma = 1.6. For our best-fitting model pa- 
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rameters, only eccentricities between 0.58 and 0.65 yield 
a mass of A between 0.5 and 1.0 Mq. Therefore, for 
a mass ratio of 1.6, the eccentricity reported by W04 
(e = 0.7) is larger than the value allowed by our radial 
velocity measurements assuming 0.5 < Ma < 1.0 Mq. 
However, this discrepancy is not too surprising since the 
W04 model used only two radial velocity measurements. 
It is also important to note that the mass limits of star A 
are derived by placing KH 15D on a theoretical HR dia- 
gram and are therefore subject to the accuracy of the 
PMS evolutionary model employed. In all other fea- 
tures of the KH 15D binary system, there is a remarkable 
agreement between the W04 model and the orbital solu- 
tion calculated from the radial velocities. Because of the 
strong evidence of a second star from our radial veloc- 
ity measurements and the historical photometry, we find 
the two-star class of model to be the most compelling ex- 
planation of the KH 15D photometric variability mecha- 
nism. 

However, the case of the "winking star" is still far from 
closed. A key missing aspect of the two-star models is 
direct detection of a circumbinary disk around KH 15D. 
He02 report a lack of near-IR excess and a null detection 
at millimeter wavelengths. CM04 state that such findings 
are consistent with a circumbinary ring having an inner 
radius of ~ 1 AU that is tidally truncated by the central 
binary, and an outer radius of ~ 5 AU that is possibly 
shepherded by an as yet unseen planet. They predict 
mid-infrared fluxes that are observable with the Spitzer 
Space Telescope. Clearly such observations will be vital 
in further development of models of the KH 15D system. 
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